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■ Abstract 

We generalize for the case of arbitrary hydrodynamical matter the 
?—i . quasi-isotropic solution of Einstein equations near cosmological sin- 

gularity, found by Lifshitz and Khalatnikov in 1960 for the case of 
' the radiation-dominated Universe. It is shown that this solution al- 

^ ■ ways exists, but dependence of terms in the quasi-isotropic expansion 

acquires a more complicated form. 

PACS: 04.20-q; 04.20.Dw 

In the paper [I]] by Lifshitz and Khalatnikov, the quasi-isotropic solution 
of the Einstein equations near a cosmological singularity was found provided 
the Universe was filled by radiation with the equation of state p — |. The 
metric of this solution was written down in the synchronous system of refer- 
ence 

ds 2 = dt 2 — r ) a pdx a dd i \ (1) 
where spatial metric 7^/3 near the singularity has the form 

7a/3 = ta a/3 + t 2 b a/3 H (2) 
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where a a p and b a p are functions of spatial coordinates. The functions a a/ 3 are 
chosen arbitrary, and then the functions b a p and also the energy and velocity 
distributions for matter can be expressed through these functions (for details 
see and also g, |). 

In correspondence with the standard cosmological model of the hot Uni- 
verse, it was supposed that the natural equation of state for the matter 
near the cosmological singularity is that of radiation: p = e/3. However, 
nowadays the situation has changed in connection with the development of 
inflationary cosmological models, which as an important ingredient contain 
inflaton scalar field or/and other exotic types of matter M. One can add 
also that the appearance of brane and M theory cosmological models || and 
the discovery of the cosmic acceleration |Bj suggests that the matter playing 
essential role on different stages of cosmological evolution can obey very dif- 
ferent equations of state M. Thus, generalization of the old quasi-isotropic 
solution of the Einstein equations near the cosmological singularity can be 
useful in this new context. 

In this note we make such a generalization for the equation of state: 

p = ke, (3) 

where p denotes pressure and e denotes energy density []. The Friedmann 
isotropic solution near the singularity for such a matter behaves as 

a ~ a t m , (4) 

where 

4 

m = — — . (5) 

3(1 + k) 1 ; 

We look for an expression for a spatial metric in the following form: 

la p = t m a aP + t n b a p, (6) 

where the power index m is given by Eq. @. We leave the power index n 
free for some time, requiring only that 

n > m. (7) 



1 Actually, the authors have known this generalization for a long time but have never 
published it in detail in regular journals. 
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The inverse metric reads 

a/3 uaf3 

^ = - (8) 

where a af3 is defined by the relation 

a a % 7 = <J« (9) 

while the indices of all the other matrices are lowered and raised by a a p and 
a a/3 , for example, 

b% = a a ^b lP . (10) 

Let us write down also expressions for the extrinsic curvature, its contractions 
and its derivatives: 

Kaf3 = ^ = mt m - l a a p + nt n ~ l b a p, (11) 
a _ m5 a , (n - m)6g 

K ° - + t m-n+l ' y LA > 

a 3m (n — m)b 
K a = —+ tm - n+l > I 13 ) 



<9fi^ m<5f (m — n + l)(n — m)6^ 
<9«;^ 3m (m — n + l)(n — m)b 

~df ~ ~~t 2 t m-n+2 ' 



(14) 
(15) 



K " K /3 - —jT + t m-n+2 " \ lKi ) 

We need also an explicit expression for the determinant of the spatial metric: 
7 = det7 a/3 = t 3m (l + r- m &)deta, (17) 

7 = ^ = (3mi 3m_1 + 6(2m + n)^"" 1 ) det a, (18) 
7 3m / b(n — m)t n ~ m \ 

7~~V + 3m" J' ' ^ 

Now, using well-known expressions for the components of the Ricci tensor 



J£=o«*-4J> (2i) 



j>P - -pP _ I^i _ SL K P (22) 

where P@ is a three-dimensional part of the Ricci tensor, and substituting 
into Eqs. ©-(|22]) the expressions ([11])- ([19]), one get 

_ 3m(2 - m) _ (n - l)(n - m)b . , 

n o - 4t 2 2t m - n+2 ' 

< = ^T(C-M, (24) 
R P = -El + m (2 " 3m 



(n — m)(2 — 2n — m)6^ m(ra — m)b5^ 



(25) 



Notice, that in Eq. ([25]) denotes a three-dimensional Ricci tensor 
constructed by using the metrics a a p. The terms in the curvature tensor P@, 
which are proportional to (3 a p have the time dependence ~ f2 J„ n and are less 
divergent than the first term in the right-hand side of Eq. (|25|) provided the 
condition ([?]) is satisfied. 

Now, let us write down the expressions for the components of the energy- 
momentum tensor of the perfect fluid 

T ik = (s + p)UiU k - pg ik , (26) 

satisfying the equation of state (H). Up to higher-order corrections, they have 
the following form : 

T ° = e (27) 

T° = e(k + IK, (28) 

TS = ~ks5i (29) 

T = T\ = e(l-3k). (30) 
Using the Einstein equations 

R{ = 8nG(Ti - ls(T), (31) 



one has from 00-component of these equations: 

1 / 3m(2 - m) (n - l)(n - m)b \ , , 

8nG£ = 3kTi{-^ — J' (32) 

and from Oa-component of these equations one has 

u _ (n-m)(3k + l)(b^-h a )t^ m 

Ua 3to(2 -m)(k + 1) ' [ ' 

Now, writing down the spatial components of the Einstein equations, using 
the expressions (|29|)-(|30D and the expression (|32] ) for the energy density e, 
one get: 

Pg m(2 - 3m)5f (n - to) (2 - 2n - m)6^ 



£m ^2 At m ~ n + 2 

m(n — m)b5^ (k — 1)5^ ( 3m(2 — to.) (n — l)(n — m)& N 
4t m ~ n+2 ~~ 3fc + 1 V t m " n + 2 



.(34) 



Using the relation (|5|) it is easy to check that the terms proportional to 
£2 in the left- and right-hand sides of Eq. ( |34"D cancel each other. On the 
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other hand, the only way to cancel the term ^ is to require that the terms 
proportional to fm _ 1 n+2 behave as the term i.e. 

n = 2. (35) 

In this case, the condition of cancellation of terms proportional to gives 
the following expression for the tensor b@: 

p AP? P^(-3m 2 + 12m-4) 

a m 2 -4 3m(m-3)(m 2 -4) ' 1 J 



Using the relation (|5]) one can rewrite the expression (|36| ) in the following 
form: 

9 ( fc + 1 ) 2 (p | (3fc 2 -6fc-5)^P \ 
(3fc + 5)(3A;+ 1) I a 9k + 5 )' { ' 



It is easy to see that the Eq. (p7\ j expressing the second-order correction 
for the spatial metric (|2]) is well defined for all the types of hydrodynamical 
matter with < k < 1, including the stiff matter, i.e. a fluid with the 
equation of state p = e. 
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Now, using the relation 



= \h, (38) 

and the formulae ([37]) and (|33|) we arrive to the following expression for the 
three-dimensional velocity u a : 

21k(k + lfP. a 3 ^_ , N 

8(3A; + 5)(9A; + 5) V ; 

Note that the velocity flow is potential. Actually, it can be shown that this 
important property remains in all higher orders of perturbative expansion for 
the quasi-isotropic solution. Similarly, the expression for the energy density 



of matter (32) can be rewritten as 



8lG£ = + »(*+D*, . (40) 

3(A; + l) 2 t 2 2 (9k + 5)t^TT) 

It is straightforward to check that further terms of perturbative expansion 
(0) for the metric j a p have the following form: 

lap = t m a af3 + t 2 b af3 + t 2+ ^c aP + ■■■. (41) 

Thus, we have seen that this expansion has a curious feature. The order of 
its first term t m is defined by the equation of state of the matter (|3|), the 
second order term always has the behavior ~ t 2 , while logarithmic distance 
between orders is equal to 2 — m. 

It is easy to understand that the quasi-isotropic expansion does work if 
and only if the first term in the right-hand side of Eq. (|6|) is smaller than 
the second one. Remembering that n = 2 for any value of m and using the 
equation (|5]), we get the following restriction on the parameter k from the 
equation of state @: 

k > ~. (42) 

Thus, for the values k < — 1/3 the quasi-isotropic expansion at small times 
may not be constructed. It is interesting to notice, however, that for k = 
const < — |, a quasi- isotropic-like solution arises as a late-time (t — > 00) 
attractor for generic inhomogeneous evolution of space-time || (this regime 
is the power-law inflation actually). Of course, perturbative expansion is 
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made in inverse powers of t in that case. Also, we would like to note that 
different aspects of relation between the quasi-isotropic expansion and other 
approximation schemes were considered in detail in the papers [ |T0|j . 

This work was supported by RFBR via grants No 02-02-16817 and 00- 
15-96699. A.K. is grateful to the CARIPLO Scientific Foundation for the 
financial support. 



References 

[1] Lifshitz E M and Khalatnikov I M 1960 ZhETF 39 149 

[2] Lifshitz E M and Khalatnikov I M 1964 Sov. Phys. Uspekhi 6 495. 

[3] Landau L D and Lifshitz E M 1979 The Classical Theory of Fields (Perg- 
amon Press) 

[4] Starobinsky A A 1979 JETP Lett 30 682; Starobinsky A A 1980 Phys. 
Lett. 91B 99; Guth A H 1981 Phys. Rev. D 23 347; Linde A D 1982 Phys. 
Lett. 108B 389; Linde A D 1983 Phys. Lett. 129B 177; Starobinsky A A 
1986 Stochastic De Sitter (inflationary) Stage in the Early Universe, in 
Field Theory, Quantum Gravity and Strings, (Eds. H.J. De Vega and N. 
Sanchez, Springer- Verlag, Berlin) 107; Linde A D 1990 Particle Physics 
and Inflationary Cosmology (Harward Academic Publishers, New York) 

[5] Banks T and Fischler W 2001 M theory observables for cosmolog- 
ical space-times, [hep-th/0102077| ; An Holographic Cosmology, |hep 

l th/onim 



[6] Perlmutter S J et al 1999 Astroph. J. 517 565; Riess A et al 1998 Astron. 
J. 116 1009 

[7] Sahni V and Starobinsky A A 2000 Int. J. Mod. Phys. D 9 373 

[8] Miiller V, Schmidt H-J and Starobinsky A A 1990 Class. Quantum Grav. 
7 1163 

[9] Lucchin F and Matarrese S 1985 Phys. Rev. D 32 1316 

[10] Comer G L, Deruelle N, Langlois D and Perry J 1994 Phys. Rev. D 49 
2759; Deruelle N and Langlois D 1995 Phys. Rev. D 52 2007 



7 



